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bo ABSTRACT 

A' 

■ Aims. We aim to use microlensing taking place in the lensed quasar QSO 2237+0305 to study the structure of the broad line region and 
| measure the size of the region emitting the C IV and C III] lines. 

^» . Methods. Based on 39 spectrophotometric monitoring data points obtained between Oct. 2004 and Dec. 2007, we derived lightcurves for 
, — . j the CIV and CIII] emission lines. We used three different techniques to analyse the microlensing signal. Different components of the lines 
[~~«, (narrow, broad, and very broad) were identified and studied. We built a library of the simulated microlensing lightcurves that reproduce the 
OO ' signal observed in the continuum and in the lines provided only the source size is changed. A Bayesian analysis scheme is then developed to 
s ! ' derive the size of the various components of the BLR. 

^sj ! Results. 1. The half-light radius of the region emitting the CIV line is found to be R CTV ~ 66+^° light-days = 0.06^^ pc = 1.7+f-J 10 17 cm 
t—{ • (at 68.3% CI). Similar values are obtained for C III]. Relative sizes of the carbon-line and V-band continuum emitting-regions are also derived 

[ with median values of j^' uc /^ cout m the range 4 to 29, depending on the FWHM of the line component. 
y—i ■ 2. The size of the C IV emitting region agrees with the radius-luminosity relationship derived from reverberation mapping. Using the virial 
1 theorem, we derive the mass of the black hole in QSO 2237+0305 to be M BH ~ ltf 3±03 M Q . 
3. We find that the CIV and CIII] lines are produced in at least 2 spatially distinct regions, the most compact one giving rise to the broadest 
component of the line. The broad and narrow line profiles are slightly different for C IV and C III]. 
<-h ' 4. Our analysis suggests a different structure for the C IV and Fe II+III emitting regions, with the latter produced in the inner part of the BLR 
or in a less extended emitting region than C IV. 



, — \ 

X 



C3 



Key words. Gravitational lensing: micro, strong, quasars: general, quasars: emission lines, quasars: individual QSO 2237+0305, line: formation 



1. Introduction 

We know that quasars and active galactic nucleii (AGN) are 
powered by matter accreted onto a supermassive black hole. 
The accretion of material in the direct vicinity of the central 



* Based on observations made with the ESO-VLT Unit Telescope 
# 2 Kueyen (Cerro Paranal, Chile; Proposals 073.B-0243(A&B), 
074.B-0270(A), 075.B-0350(A), 076.B-0197(A), 177.B-0615(A&B), 
PI: F. Courbin). 



black hole releases most of the quasar energy in the form of 
power-law continuum emission. Ionised gas surrounds the cen- 
tral accretion disc and gives rise to broad emission lines, which 
are used as footprints that allow the identification and classifi- 
cation of quasars. Our knowledge of the kinematics and phys- 
ical conditions prevailing in the BLR gas remain elusive, es- 
pecially because the nuclear region of quasars is still spatially 
unresolved with existing instrumentation. 
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Current insights into the BLR come from various kind of 
studies: empirical modelling of the line shape with kinemati- 
cal models, use of photo-ionisation codes to reproduce the ob- 
served flux ratios between spectral lines, spectropolarimetric 
observations, statistical study of the width and asymmetry of 
the lines, use of the principal component analysis technique, 
and velocity resolved reverberation mapping (e.g. Boroson & 
Green [19921 Sulentic et al. [20001 Smith et al. [2005] Marziani 
et al. 120061 Zamfir et al. [20081 Gaskell [20091 1201 Obi Bentz 
et al. 1201 01) . Despite the development and many successes 
of these methods, as briefly summarised below, we still do 
not completely understand the structure and kinematics of the 
BLR. The microlensing of the broad emission lines in multi- 
ply imaged lensed AGNs provides us with a powerful alterna- 
tive technique for looking at the BLR, measure its size even in 
high luminosity distant quasars, and get hints of the structure 
and geometry of both emission and intrinsic absorption within 
the BLR (e.g. Schneider & Wambsganss 1990 Hutsemekers 
et al. [19941 Lewis & Belle [T9981 Abajas et al. [20021 Popovic 
et al. 120031 Lewis & Ibata 120041 [20061 Richards et al. [20041 
Abajas et al. 120071 Sluse et al. [20071 [20081 Hutsemekers et 
al. l2%T0l) . 



1.1. Phenomenology of the line profiles 

Our primary clue on the BLR comes from the shape of the 
broad emission lines. The most detailed studies of broad emis- 
sion lines have focused on two lines: H/?/l4863 and CIV /11549 
(Sulentic et al. I2000I for a review). Of direct interest for the 
present work is C IV/11549. The C IV profile shows a broad va- 
riety of shapes from strongly asymmetric to symmetric (e.g. 
Wills et al. [19931 Baskin & Laor I20051 has an equivalent 
width anti-correlated with its intensity (Francis et al. 1992 ), and 
shows greater variability in the wings than in the core (Wilhite 
et al. 2006 ). In addition, the C IV line is also found to be sys- 
tematically blueshifted by several hundred to a few thousand 
kilometres per second compared to the low ionisation lines 
(Gaskell [191321 Corbin [19901 Vanden Berk et al. [2000T >. The 
analysis of about 4000 SDSS quasars by Richards et al. (|2002) 
suggests that this shift is caused by a lack of flux in the red wing 
of the line profile and correlates with the quasar orientation. 
Variability studies and emission-line decomposition techniques 
(principal component analysis and analytical fitting) indicate 
that the region emitting C IV could be built up with two com- 
ponents, a "narrow" emission core of FWHM ~ 2000 kms 1 
emitted in an intermediate line region (ILR) possibly corre- 
sponding to the inner part of the narrow line region and a very 
broad component (VBC) with FWHM -7000 kms" 1 produc- 
ing the line wings (Wills et al. 119931 Brotherton et al. 1994a 
Sulentic 2000, Wilhite et al. I2006I I. In radio-quiet objects, the 
VBC is observed to be systematically blueshifted by thou- 
sands of kms 1 with respect to the narrow core (Brotherton 
et al. I1994al Corbin [1995) suggesting it is associated with 
outflowing material. The ILR component disclosed in C IV is 
probably different from the one recently uncovered in H/3 (Hu 
et al. 120081 1. but its exact nature and the physical conditions 



in this region are still being debated (Brotherton et al. 1994b 
Sulentic & Marziani [19991 Marziani et ai. l20"0"6l 

The properties of the C III] emission line have received less 
attention in the literature, probably mostly because the line is 
blended with A1III Al 857, Si III] ^1892 and, anFeII+IIU1914 
complex. In their study of the CIV and CIII] emission, 
Brotherton et al. (|1994al > find that the 2 lines often have dif- 
ferent profiles and that, for some objects, a two-component de- 
composition (i.e. narrow core+very broad wings) does not pro- 
vide a good model of C III], so a third component is required. 
They also find that the VBC needed to reproduce C III] has to 
be larger than the corresponding component of C IV. This sup- 
ports the idea that the region emitting C III] is different from 
the one emitting CIV (Snedden & Gaskell [T^99l l. 

Recently, Marziani et al. (120101 1. based on the line profile 
decomposition of a small sample of AGNs selected in the 4D 
Eigenvector 1 context (e.g. Boroson & Green [19921 Zamfir et 
al. [2008 ), suggest that all the broad emission line profiles are 
composed of three components of variable relative intensity 
and centroid shift (from line to line in a given object and be- 
tween objects). They suggest a classical unshifted broad com- 
ponent (FWHM=600-5000 krns -1 ), a redshifted very broad 
component and a blueshifted component mostly visible in the 
so-called Population A objects (i.e. objects with FWHM < 
4000 kms~'). Based on line ratios, they also tentatively infer 
that these components arise from different emitting regions. 

The smoothness of the line profiles and the physical condi- 
tions derived from photo-ionisation models allowed several au- 
thors to put constraints on the "structure" of the BLR gas. Two 
popular models remain. The first one considers that the BLR is 
a clumpy flow composed of small gas clouds, and the second 
one assumes a smooth gas outflow originating in an accretion 
disc (Elvis 120001 Laor I2007l and ref. therein). Several geome- 
tries have been considered for the BLR gas, the most popular 
ones being disc-like, spherical and biconical models (e.g. Chen 
& Halpern [19891 Robinson [T995l l. Murray and Chiang ( fT997T > 
demonstrate that a continuous wind of gas originating in an ac- 
cretion disc can successfully reproduce the profile and system- 
atic blueshift of the C IV emission line. Other indications that 
a fraction of the BLR material has a disc-like geometry comes 
from statistical studies of the broad line profiles in samples of 
(mostly radio-loud) AGNs (Vestergard et al. 120001 McLure 
& Dunlop [20021 Jarvis & McLure [20061 Decarli et al. [20081 
Risaliti et al. 2010), and from spectro-polarisation observations 
of Balmer lines in AGNs (Smith et al. 120051 1. Some of these 
studies also suggest there is a second, spherically symmetric 
component with Keplerian motion. The spectropolarimetric ob- 
servations of PG 1700+518 by Young et al. (120071) support a 
disc+wind model for the Hff emission. On the other hand, bi- 
conical models of the BLR seem needed to explain the vari- 
ability of the rare double-peaked AGNs (i.e. AGNs showing 
emission lines with two peaks) and their polarisation properties 
(Sulentic et al. 119951 Corbett et al. 119981 ). These studies show 
that there is no consensus on the geometry of the BLR. From an 
observational and theoretical perspective, it seems nevertheless 
to depend on the ionisation degree of the line and on the radio 
properties of the object. 
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1.2. The radius-luminosity relationship 

The "size" of the the broad line region Rblr has been measured 
in about 40 AGNs by use of the reverberation mapping tech- 
nique (e.g. Rrolik et al. I1991I Home et al. 2004). The empir- 
ical relation Rblr k L a (a -0.5-0.6, Kaspi et al. l2005l Bentz 
et al. 12006b . combined with the virial theorem, allows one to 
derive a relation linking the black hole mass Mbh, the AGN lu- 
minosity L and the FWHM of the emission line. This relation 
is one of the most popular methods used to measure black hole 
masses based on single epoch spectroscopic data and study 
their growth, evolution and correlation to other AGN proper- 
ties. 

The Rblr 00 L a relation has been derived quite accurately 
for the broad component of the H/3 line, but only a few ob- 
jects have Rblr measurements for high ionisation lines like 
CIV (e.g. Kaspi et al. 120071 1. This is a severe problem for 
black hole mass measurements of high redshift objects. The 
use of the C IV line to derive black hole masses is desirable 
but faces several problems related to our understanding of the 
structure and geometry of the region emitting this line (see e.g. 
Marziani et al. 120061 for a review). The existence and con- 
tamination of a narrow emission component in CIV which 
could bias FWHMc iv (Bachev et al.[2004 ) and the possible ab- 
sence of virial equilibrium, especially in sources showing large 
blueshifts (Richards et al. 120021) . are major concerns. 

Velocity resolved reverberation mapping (Home et 
al. [2004 ) should allow one to get insights on these problems, 
but the technique is still under development and has not yet 
been applied to the C IV line. Most of the recent advances in 
velocity resolved echo-mapping are based on the analysis of the 
H/3 emission line (Denney et al. |2009l Bentz et al. |20U9ll20T0l l. 
Although the kinematics of the Balmer gas was thought to be 
relatively simple, the technique provided puzzling results as 
both keplerian rotation, inflow and outflow signatures appear 
in various objects. It is still unclear whether this reveals a wide 
variety of kinematics, a biased interpretation of the observed 
signal, or a spurious signal introduced by observational arti- 
facts. The past evidence of inflow and outflow signatures in 
NGC 5548 (Clavel et al. fT99T1 Peterson et al. fl99Tl Kollatschny 
& Dietrich 1996 ) might indicate that inflow/outflow signatures 
are not unambiguous. A possible explanation of the variety 
of observed signals might be off-axis illumination of the BLR 
(Gaskell 2010a 2010b and references therein). 

1.3. Microlensing in QSO 2237+0305 

The many open questions concerning the BLR we outlined 
above motivate the interest in developing new techniques of 
probing this region. In this paper we investigate the con- 
straints on the broad lines provided by means of the mi- 
crolensing study of the lensed quasar QSO 2237+0305, for 
which long-term spectro-photometric monitoring has been car- 
ried out (Eigenbrod et al. l2007l hereafter Paper /). The gravita- 
tional lens QSO 2237+0305, also known as "Huchra's lens" 
or the "Einstein Cross", was discovered by Huchra et al. 
(1985) during the Center for Astrophysics Redshift Survey. 
It consists of a z s = 1.695 quasar gravitationally lensed into 



four images arranged in a cross-like pattern around the nu- 
cleus of a zi = 0.0394 barred Sab galaxy. The average pro- 
jected distance of the images from the lens centre is only 
700 pc ~ 0.9", such that the matter along the line of sight to 
the lensed images is mostly composed of stars. The symmet- 
ric configuration of the lensed images around the lens-galaxy 
bulge ensures a time delay of less than a day between the 
lensed images, such that intrinsic flux variations should be 
seen quasi- simultaneously in the four images (Rix et al. 1992, 
Wambsganss & Paczynski 1994). Since the low redshift of 
the lens galaxy leads to a high relative transverse velocity be- 
tween the observer, the lens, and the source, the lensed im- 
ages of QSO 2237+0305 are continuously flickering due to mi- 
crolensing produced by the stars in the lens galaxy, on smaller 
timescales than in any other lens. The microlensing affecting 
the images of QSO 2237+0305 leads to large variations in 
amplitude (Udalski et al. 120061) . reaching up > 1 mag, often 
accompanied by chromatic microlensing of the quasar contin- 
uum (Wambsganss & Paczy nski l 1 99 fb used to study the accre- 
tion disc temperature profile (e.g. Kochanek 2004, Anguita et 
al.[2p08, Eigenbrod et al. 2QQ% -Paper IT). Several studies have 
shown that the broad emission lines are also significantly af- 
fected by microlensing (Metcalf et al. 120041 Wayth et al. 120051 
Paper I). Based on the microlensing of CHI] observed at one 
epoch in 2002, Wayth et al. (2005) derived a most likely size of 
the region emitting this line of 0.06 h 1 ^ 2 pc. The work presented 
here aims at improving the measurement of the BLR size by us- 
ing monitoring data instead of a single epoch measurement and 
at constraining the structure of the BLR. 

The structure of the paper is the following. In Sect. [2] we 
briefly summarise the data we used and explain the three dif- 
ferent methods we applied to analyse the spectra. In Sect. [3j 
we apply these techniques to our data to understand how mi- 
crolensing deforms the emission lines and to measure the flux 
ratios in several portions of the C IV and C III] emission lines. 
In Sect. [4] we explain the microlensing simulations we devel- 
oped to derive the size of the broad emission lines and contin- 
uum emission region. In Sect. [5] we present our results: mea- 
surement of the size of the BLR and of the continuum region, 
comparison of our BLR size with reverberation mapping, esti- 
mate of the black hole mass, and finally constraints on the BLR 
structure. Finally, we summarise our main findings in Sect. [6] 

2. Observations and analysis techniques 

We compare emission line fluxes in images A & D of 
QSO 2237+0305. A comprehensive description of the obser- 
vations and of the observational setup can be found in Papers I 
& II. The full data set consists of a series of spectra obtained at 
39 different epochs between October 2004 and December 2007 
with the FORS 1 instrument at the ESO Very Large Telescope. 
All the observations were carried out in MOS (Multi Object 
Spectroscopy) mode. The multi-object mask used to obtain the 
data presented here included the quasar images A and D, as 
well as several field stars which are used for simultaneous flux 
calibration and deconvolution. The spectra of the individual 
lensed quasar images were de-blended (from each other and 
from the lensing galaxy) using the MCS deconvolution tech- 
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Fig. 1. Comparison of V-band OGLE lightcurves of images A 
& D (points) of QSO 2237+0305 with the lightcurves derived 
from our spectra (diamonds; crosses are used for epochs devi- 
ating from the OGLE lightcurves). Error bars (not displayed) 
are typically of the order of 0.02 mag (A) and 0.04 mag (D). 



nique (Magain et al. 119981 ). A slit width of 0.6"was used to- 
gether with grism GRIS_300V and blocking filter GG375, lead- 
ing to a resolving power R = 400 at 5900 A. For the data ob- 
tained in 2007, the FORS CCD had been replaced by a new 
camera that was more sensitive in the blue wavelength range. 
With this new camera, the GG375+30 blocking filter was not 
used. The spectral resolution was improved by about 10%, but 
the observations obtained with the new setup have a signal-to- 
noise ratio S/N<10 in the red wavelength range {A > 6700 A) 
due to fringing. Table Q] provides a log of the data sample we 
used. Data flagged with a f are systematically removed as the 
C IV emission is not covered by our spectra. 

We show in Fig.Q]a comparison of the V-band lightcurves 
of images A & D provided by the OGLE collaboration (Udalski 
et al. 120061 1 with the lightcurves derived by applying V-band 
response curve to the spectro-photometric monitoring data. 
Throughout this paper, we shift the magnitude of image D from 
the OGLE-III lightcurve by -0.5 mag with respect to the pub- 
lished values to account for different flux calibration of that im- 
age in the OGLE-III data (see Paper I for discussion). Although 
the overall agreement is very good, several epochs identified 
with an 'X' symbol deviate from the OGLE lightcurve. These 
data are flagged with an asterisk (*) in Tableland we test how 
they affect our results in Sect. 15.21 

To study how microlensing deforms our spectra and how 
the broad lines are affected, we use three different tech- 
niques that are explained hereafter. First, we apply the multi- 
component decomposition technique (MCD); i.e., we assume 
that multiple components cause the quasar emission and mea- 
sure the flux in these individual components (Sect. 12. U . 
Second, we use the macro-micro decomposition technique 
(MmD); i.e., we follow the prescription of Sluse et al. ([2007) 
and Hutsemekers et al. (2010) to decompose the quasar spec- 



Table 1. Journal of the observations used in this paper. The 
observation ID refers to the whole data set presented in Papers 
I and II. Data flagged with a t are not used. Data flagged with 
an * are those deviating from the OGLE lightcurve in Fig. [1] 
As in Paper II, epochs #7, 8, 10, 30 have been excluded from 
the analysis. 





ID 


Civil Date 


JD-2450000 


Seeing ["] 


Airmass 




1* 


2004-10-13 


3292 


0.86 


1.204 




2 


2004-11-14 


3324 


0.75 


1.184 




3 


2004-12-01 


3341 


0.88 


1.355 




4* 


2004-12-15 


3355 


0.99 


1.712 




5 


2005-05-11 


3502 


0.87 


1.568 




6 


2005-06-01 


3523 


0.63 


1.342 


IX 


9 


2005-08-06 


3589 


0.51 


1.135 


73 
O 


11 


2005-08-25 


3608 


0.49 


1.261 


■c 

Q 


12 


2005-09-12 


3626 


0.70 


1.535 


0- 


13* 


2005-09-27 


3641 


0.92 


1.480 




14* 


2005-10-01 


3645 


0.78 


1.281 




15 


2005-10-11 


3655 


0.57 


1.140 




16* 


2005-10-21 


3665 


0.70 


1.215 




17* 


2005-11-11 


3686 


0.90 


1.137 




18* 


2005-11-24 


3699 


0.78 


1.265 




19* 


2005-12-06 


3711 


1.10 


1.720 




20+ 


2006-05-24 


3880 


0.87 


1.709 




21+ 


2006-06-16 


3903 


0.66 


1.213 




22 1 "* 


2006-06-20 


3907 


0.64 


1.286 


<N 
Oh 


23 1 " 


2006-06-27 


3914 


0.41 


1.145 


24 


2006-07-27 


3944 


0.74 


1.316 


Perioc 


25' 


2006-08-03 


3951 


0.73 


1.246 


26 


2006-10-13 


4022 


0.59 


1.176 




27 


2006-10-28 


4037 


0.57 


1.148 




28 


2006-11-10 


4050 


0.89 


1.515 




29 


2006-11-27 


4067 


0.87 


1.255 




31 


2006-12-22 


4092 


0.80 


2.018 




32* 


2007-07-10 


4292 


0.63 


1.153 




33* 


2007-07-15 


4297 


0.57 


1.158 




34 


2007-07-25 


4307 


0.79 


1.231 




35 


2007-08-03 


4316 


0.68 


1.412 


0- 


36 


2007-08-27 


4340 


0.88 


1.133 


o 


37* 


2007-09-06 


4350 


0.67 


1.396 


•e 

u 


38 


2007-09-20 


4364 


0.73 


1.230 


Oh 


39* 


2007-09-23 


4367 


1.23 


1.530 




40 


2007-10-05 


4379 


0.59 


1.153 




41 


2007-10-10 


4384 


0.76 


1.283 




42 


2007-11-15 


4420 


1.07 


1.189 




43* 


2007-12-01 


4436 


1.00 


1.502 



trum into a component Fm, unaffected by microlensing and an- 
other one, Fm/j, which is microlensed (Sect. 12.21 ). Third, we 
use the narrow band technique (NBD); i.e., we cut the emis- 
sion lines into different velocity slices and measure the flux 
in each slice to look for microlensing within the emission line 
(Sect. 1231). 
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2.1. The multi-component decomposition (MCD) 

Our most comprehensive analysis technique is a multi- 
component spectral decomposition of the quasar spectrum 
(MCD). With this technique, we decompose the quasar spectra 
into the components likely associated with different emission 
regions (Dietrich et al. 2003 ). As in Papers I & II, the spectrum 
is modelled as a sum of three components: i) a power-law for 
the quasar continuum emission, ii) an empirical iron template, 
and iii) a sum of Gaussian profiles to model the emission lines. 
The exact decomposition of each spectrum differs in several 
respects from the one presented in Papers I & II: 

- First, the quasar emission in image D is corrected for dust 
reddening produced by the lensing galaxy using the dif- 
ferential extinction values derived in Paper I. We adopt a 
Cardelli (1989) Milky-way like extinction law with (Rv, 
Av) = (3.1, 0.2 mag). 

- Because the UV Fell+III template constructed by 
Vestergaard & Wilkes (2001) from I Zwl does not per- 
fectly reproduce the Fell+III blended emission observed 
in QSO 2237+0305 , we construct the iron template em- 
pirically. In order to take possible intrinsic flux variations 
and/or systematic errors into account, we proceed indepen- 
dently for data obtained in 2004-2005 (period PI), 2006 
(period P2), and 2007 (period P3). This template is con- 
structed iteratively based on the residuals of the fitting pro- 
cedure. In a first step, we fit only a power-law to the spec- 
tra, excluding those regions where Fe II+III is known to be 
significant (Vestergaard & Wilkes 120011 ) and where there 
are noticeable quasar broad emission lines. The subsequent 
fits are performed on the whole quasar spectrum (exclud- 
ing atmospheric absorption windows), fixing the width and 
centroid of the quasar emission lines. A new fit is then per- 
formed using the residuals as the iron template. The proce- 
dure is repeated until a good fit of the whole quasar spec- 
trum is obtained. 

- In order to better reproduce the detailed shape of the C IV 
and CIII] line profiles, we refine their modelling com- 
pared to Paper I by decomposing each line profile into a 
sum of three Gaussian profiles. The C IV line is composed 
of two emission (a narrow-component NC and a broad- 
component BC) and of one absorption component (AC). 
The faint 1600 A blend of emission sometimes observed 
in AGNs, and caused mostly by HeIU1640, O III] ,11 663 
(Fine et al. 120 101 and ref. therein), is included by construc- 
tion in our Fe II+III template. To reproduce the C III] line, 
we need to use three emission components of increasing 
width (abbreviated NC, BC, and VBC). With the exception 
of CIII] (VBC), the central wavelength A c of the Gaussians 
are allowed to vary from epoch to epoch as is the ampli- 
tude of each component. Table |2]provides for each image, 
the central wavelength A c and the FWHM of the different 
components of the line. The average value and standard 
deviation between epochs are displayed. When the stan- 
dard deviation is not displayed, it means that the quantity 
is fixed for all epochs (see Appendix lAl for details). Figure 
|2] shows the CIV and CIII] emission line on 2004-11-14 
and the corresponding profile decomposition. As in Paper 
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Table 2. Decomposition of the C IV and C III] emission lines 
(Average value and scatter between epochs when the parameter 
is fitted). The central wavelength A c is expressed in A and the 
FWHM inkms -1 . 



ID A D 

I c FWHM T c FWHM 
CIV 

~BC 1547.8 ± 1.0 6332 ± 404 1546.9 ± 0.9 6128 + 557 

AC 1548.9 ±0.5 1 153 ± 190 1548.7 ±0.4 1 186 ± 253 

NC 1547.6 ± 0.6 2580 1547.0 ± 0.4 2580 

CIII] 

BC 1910.7 ± 0.6 3400 1910.7 ± 1.1 3400 

NC 1907.3 ± 0.4 1545 1907.5 ± 0.3 1545 

VBC 1910.9 8548 1910.9 8548 



I, the other emission line profiles (e.g. Al III, Hell, Si III], 
OIII], Mgll) are fitted with single Gaussian profiles of 
fixed width. 

The results of the MCD decomposition are shown and 
discussed in Sect. 13.11 One could suspect that residual flux 
from the lensing galaxy is included in our empirical Fell tem- 
plate. We tested this possibility by fitting a lens galaxy tem- 
plate (as retrieved from the deconvolution of the 2D spec- 
trum, see Paper I) as an independent component. However, the 
quasar lightcurves derived in this case strongly disagree with 
the OGLE lightcurves, indicating that residual contamination 
of the spectra by the lensing galaxy is not a big issue. Finally, 
we emphasise that the decomposition of the carbon lines in 
multiple Gaussians of different widths is empirical and does not 
necessarily isolate emission lines coming from physically dif- 
ferent regions. The adopted naming of the components is based 
on their width, and is line dependent and internal to this paper. 
Because of the different decompositions found for the C IV and 
C III] lines, comparison between Gaussian components should 
be performed with care, keeping the actual line widths in mind. 

2.2. Macro-micro decomposition (MmD) 

This technique, introduced in Sluse et al. (2007]) and 
Hutsemekers et al. (2010) and previously nicknamed Fm-Fm^ 
decomposition, allows us to determine the fraction of the 
quasar emission that is not affected by microlensing indepen- 
dently of any modelling of the quasar spectra. Because spectral 
differences between lensed images only stem from microlens- 
ing (here we correct for differential extinction prior to the de- 
composition), we can assume that the observed spectrum F, 
of image i is the superposition of an emission Fm that is only 
macro-lensed and of another component Fm/j, both macro- and 
m/cro-lensed. Using pairs of spectra of different lensed im- 
ages, it is then easy to extract both components Fm and F^. 
Defining M = Ma/Md(> 0) as the macro-magnification ratio 
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between image A and image D, and /j. = ^//io as the mi- 
crolensing ratio between image A and image D, we have 

F A =MxF M + MxnxF MfI 
Fd = F M + Fm^ ■ 

These 2 equations can be rewritten as 

F _ F A jM-iixF D 



where /j must be chosen to satisfy the positivity constraint 
Fm > and Fmh > 0. For the macro-magnification factor M, 
we adopt the ratio M = Ma /Mo = 1.0 provided by the lens 
model of the Einstein cross performed by Schmidt et al. ( 1998 ). 
The microlensing ratio fi(t, A) depends on wavelength and time 
and can be estimated for each epoch based on the ratio of the 
continuum emission between A & D (Sect. 12. U . Our choice 
of the continuum flux ratio to derive jj(t, A) means that we im- 
plicitly assume that the emission line and pseudo-continuum 
emission arise from regions equally or less microlensed than 
the continuum, in agreement with theoretical expectations. This 
assumption is also verified a posteriori as our analysis confirms 
that the continuum emission is also the most microlensed. 

2.3. Narrow band technique (NBD) 

A simple method of identifying differential microlensing 
within an emission line is to divide the line into various wave- 
length bands and calculate the flux ratio Fa/Fd for each of 
these bands. For the C IV and C III] emission, we chose three 
different wavelength ranges corresponding to the blue wing 
(BE; -8000 < v < -1500 km s" 1 ), the line core (CE; -1500 < v < 
1500 km s 1 ), and the red wing (RE; 1500 < v < 8000 km s 1 ) 
of the line. Before integrating the flux in these ranges, we 
first remove the continuum and the Fell+III emission from 
the spectra based on the results of the MCD fit (Sect. [XT). For 
the C III] emission, we also remove the Al III and Si III] emis- 
sion. The central wavelength corresponding to zero velocity is 
1549 A for C IV and 1907.5 A for C III]. 

3. Microlensing variability in the emission lines 

In Paper I, we reported evidence that the broad emission lines 
in QSO 2237+0305 are significantly microlensed in image A 
during the monitoring period. In this section, we refine and ex- 
tend these results (including 2007 data, which were absent from 
Paper I) using the techniques outlined in Sect. [2] This will al- 
low us to derive lightcurves for the C IV and C III] emission 
lines and identify signatures of microlensing. 

3.1. CIV and C III] emission lines 

First, we investigate how microlensing affects the emission line 
shape. For this purpose we use the MmD technique (Sect. l2~!2l . 
The results of this decomposition for the C IV and C III] emis- 
sion are shown in Fig. [3] for four representative epochs. In this 
figure, we also overplot the continuum model derived from the 



MCD decomposition and the spectrum of image D. As indi- 
cated by Eq.Q] the spectrum of image D should match Fj^ in 
the region of the spectra which are microlensed like the con- 
tinuum (Fm m = Fd when Fm = 0). The match between Fr> 
and Fnift in the wings of the emission lines reveals that the 
corresponding emitting region is microlensed as much as the 
continuum (Figure [3). It is also apparent that the microlensed 
fraction of the lines has a broad Gaussian shape with no flux 
depletion or flat core in the centre (the depletion observed for 
the CIV emission is caused by the intrinsic absorber). On the 
other hand, in Fm, only a narrow emission line is visible for 
C IV and C III]. If we remove the continuum power-law model 
and Fe II model prior to the decomposition, and plot the C III] 
and C IV lines on a velocity scale (Figure @}, we find that the 
decomposed line profiles are rather similar. The remainder of 
the difference could be caused by Fe II blended with the C III] 
emission, hence not included in our Fell template, but may 
also be real. It is also apparent that the line profile is asym- 
metric with a stronger blue component and that the C IV line is 
slightly blueshifted (~500kms _1 ) compared to C III]. Although 
comparison between epochs is difficult due to intrinsic variabil- 
ity of the emission lines, we observe that the decomposition 
remains roughly the same during the period PI (not shown). 
During periods P2 and P3, both fractions Fm and Fm^ of the 
carbon lines follow the intrinsic brightening of the continuum. 
We also tentatively see that, during period P3, the line wings 
appear in Fm, indicating that the broadest component is differ- 
ently microlensed than the continuum at these epochs. We note 
that the use of a macro-magnification ratio M ~0.85, closer to 
the radio and mid-infrared flux ratios (Falco et al. 1 1 9961 Agol 
et al. 2009, Minezaki et al. 2009), does not change the general 
features observed in Figure [3] 

The slicing of the emission lines in velocity also allows us 
to study differential microlensing in the emission lines indepen- 
dently of a model of the line profile. Figure [5] displays the flux 
measured in the blue wing (BE), red wing (RE), and line core 
(CE) of the carbon emission lines, using the NBD (Sect. l2~3l >. 
Standard error bars are calculated based on the photon noise. 
Although the error bars are large, we clearly see a difference 
of about 0.3 mag between the flux ratios measured in the core 
and the wings of the emission lines. Table[3]shows the average 
flux ratios for periods PI, P2, P3. Error bars for a period are 
standard errors calculated based on the dispersion of the data 
points. The systematic difference observed between the blue 
and red wings of the C IV profile seems to be mostly caused 
by an asymmetry of the line's core. Indeed, the offset between 
blue and red wing decreases as we increase the width of the 
central velocity range (not shown). The asymmetry of the Fm 
fraction of the line visible in Figure [4] supports this finding. 
When looking to the intrinsic variability of the different bands 
(i.e. individual lightcurves for images A & D), we observe that 
the wings vary more significantly on short timescales than does 
the core of the line. 

The results of the MCD decomposition are presented in 
Figure [6] They broadly agree with the velocity slicing tech- 
nique (Pig. [21 Table [3). Figure [6] shows that the flux ratio in the 
full C III] profile varies weakly over the monitoring interval. A 
similar flux ratio is measured in C IV and C III] with evidence 
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Fig. 2. Example of multi component decomposition of the C IV (left) and C III] {right) line profiles (Sect. 12.11) after continuum 
and Fe II subtraction. The narrowest component of the line (NC) is displayed with a dotted-red line, the broad component (BC) 
with a dashed-blue line, the very broad component (VBC) of C III] and the absorber in front of C IV are shown in green dot- 
dashed. The sum of the individual components, corresponding to the line model, is shown with a solid grey line. Characteristics 
of the components are provided in Table [2] 



CIV Gil] + Si III] CIV Clll] + Si III] 




\,„ (A) Ke*t (A) \„ (A) \,„ (A) 



Fig. 3. Macro-micro decomposition (MmD) of the C IV {left) and C III]+Si III] (right) emission lines obtained from the spectra 
of images A & D at 4 different epochs. In each panel, the black solid line shows the fraction of the spectrum Fm m affected by 
microlensing and the red solid line shows the emission Fm which is too large to be microlensed. For comparison, we also display 
the observed spectrum of image D with a dotted blue line and the power-law continuum used to calculate fi(A, t) with a dashed 
blue line. 



of a small increase (Aihad closer to 0) in 2007. Figure [6] also 
confirms that the flux ratio measured in the narrower compo- 
nent of the emission lines is closer to 1 (i.e. the macro-lensed 
ratio Ma/Mo), while the flux ratio amounts to nearly 2.5 in the 
broadest component of the lines (i.e. BC of C IV and VBC of 

cm]). 



3.2. Other emission lines 

An accurate study of the microlensing occurring in the 
other emission lines is very difficult because of the pseudo- 
continuum Fe II+III in emission blended with these lines. The 
MmD technique gives qualitative information about microlens- 
ing for the other emission lines. Figure|7]shows the Fm and Fm^ 
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Fig. 4. Macro-micro decomposition (MmD) technique applied to the C IV (black) and C III] (red) emission lines observed on 
2006-10-13. Left: Fraction of the flux Fm non affected by microlensing, Centre Fraction of the line Fmh affected by microlensing, 
Right: Full emission line profile. The intensity of C III] are multiplied by 2.2 to ease visual comparison between C IV and C III] 
profiles. The vertical dashed line indicates the velocity zeropoint corresponding to the systemic redshift. 




OGLE continuum 
Continuum A(C1V) 




JD-2450000 (days) 



4200 4400 4600 



Fig. 6. Time variation of the magnitude difference between images A & D, Amad, as measured in the continuum (centre) and in 
the different components (Fig. |2j of the CIV (left), CIII] (right) emission lines. The shaded area indicates the range of Amad 
estimated from the macro-model and from MIR measurements. The acronyms NC, BC and VBC refer to the different components 
(of increasing width) of the C IV & C III] lines as shown in Fig. [2] 



Table 3. Average values of magnitude difference Aihad as mea- 
sured in the blue wing (BE), line core (CE), and red wing (RE) 
of the C IV and C III] emission lines. 



Am AD (mag) 


Line 


Period 


BE 


CE 


RE 


CIV 


PI 


-0.91 ± 0.02 


-0.66 ± 0.02 


-1.27 ± 0.02 




P2 


-0.80 ± 0.03 


-0.66 ± 0.02 


-1.21 ± 0.04 




P3 


-0.66 ± 0.02 


-0.58 ± 0.02 


-1.10 ± 0.02 


CIII] 


PI 


-1.08 ±0.02 


-0.59 ± 0.01 


-1.07 ± 0.02 




P2 


-1.18 ±0.03 


-0.62 ± 0.01 


-1.14 ± 0.03 




P3 


-1.32 ± 0.05 


-0.61 ± 0.02 


-1.13 ± 0.02 



spectra averaged over period PI . We clearly see that not only do 



C IV and C III] have their broadest component microlensed, but 
also MgIU2798, A1IU1671, HeIU1640, and A1IIU1857. 
The narrow component of these lines is clearly visible in 
Fm- Interestingly, the Si III] /U892 is not microlensed. The ab- 
sence of microlensed Si III] emission is also clear in Fig. [3] 
The bottom panel shows the Fe^y template of Vestergaard & 
Wilkes d200Tb convolved with a Gaussian of FWHM ~ 2000 
kms" 1 . Although the iron emission in QSO 2237+0305 dif- 
fers slightly from the Vestergaard & Wilkes (2001) template 
(Paper I, Sect. 12. U . we clearly see that only the microlensed 
pseudo-continuum (i.e. emission above the continuum in F^) 
resembles the Feuv emission but, there is no sign of UV iron 
emission in Fm- The pseudo-continuum emission visible in Fm 
at A > 1920 A is unlikely to be produced by Feuv because its 
shape deviates significantly from expectations. The most likely 
origin(s) of this emission is (are) UV emission from the host 
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Fig. 5. Time dependence of the magnitude difference between 
images A & D, Am^D, as measured in the blue wing (BE; blue 
diamond), line core (CE; green square) and red wing (RE; red 
hexagons) of the C IV (Top) and C III] (bottom) emission lines. 
The shaded area indicates the range of Am^o estimated from 
the macro-model and from MIR measurements. 



galaxy, Balmer continuum emission, and scattered continuum 
light. Although this feature peaks at about 2200 A in the quasar 
rest frame, it is unlikely to be caused by differential extinction 
in the quasar host galaxy because the light-rays separation be- 
tween A & D in the host is too small (Falco et al. 1999). 



4. Microlensing simulations 

In order to derive the size of the regions emitting the car- 
bon lines, we must compare the observed microlensing signal 
with simulated microlensing lightcurves generated for differ- 
ent source sizes. The strategy we adopt is similar to the one 
explained in Paper II. Instead of searching for the best sim- 
ulated lightcurve reproducing the data, we follow a Bayesian 
scheme for which a probability is associated to each simulated 
lightcurve. A probability distribution is then derived for each 
parameter that describes the lightcurve. We provide a descrip- 
tion of the crucial steps in the next sections. 




1600 1800 2000 2200 2400 2600 2800 3000 

Km (A) 

Fig. 7. Top: Macro-micro decomposition (MmD) applied to 
the spectra of images A & D averaged over period PI. The 
black solid line shows the fraction of the spectrum Fm^ af- 
fected by microlensing, and the red solid line shows the emis- 
sion Fm, which is too large to be microlensed. For compari- 
son, we also display the observed spectrum of image D with 
a dotted-blue line and the power-law continuum of D with a 
dashed black line. Bottom: Fe II+III template from Vestergaard 
& Wilkes I200T1 convolved with a Gaussian of FWHM=2000 
kms -1 . 



4.1. Creating the microlensing pattern 

We use the state-of-the-art inverse ray-shooting code devel- 
oped by Wambsganss (1990, 1999, 2001) to construct micro- 
magnification maps for images A & D. To create these patterns, 
the code follows a large number of light-rays (of the order of 
10 10 ), from the observer to the source plane, through a field of 
randomly distributed stars. The surface density of stars and the 
shear y at the location of the lensed images are those provided 
by a macro model of the system. Instead of the surface den- 
sity, we used the dimensionless quantity k (the convergence), 
defined as the ratio between the surface density and the critical 
density. We use (k, y) A = (0.394, 0.395) and (*, y) D = (0.635, 
0.623) (Kochanek 2004). Due to the location of the lensed im- 
ages behind the bulge of the lens, we assume that 100% of 
the matter is in the form of compact objects. This assumption 
is motivated by the results of Kochanek (2004), who demon- 
strate that the most likely K smooth = in QSO 2237+0305. The 
mass function of microlenses has little effect on the simula- 
tions (Wambsganss [T992"l Lewis & Irwin [T9951 [19961 Congdon 
et al. 2007]), so we assume identical masses in the simulations. 
The mass of the microlenses sets the Einstein radius r^. For 
QSO 2237+0305, the Einstein radius projects onto the source 
plane as 



I'E 



1^^=9.84X10- 



M 



yj 0.3M© 



cm, 



(3) 



where the D are angular diameter distances, and the indices o, 
I, s refer to observer, lens and source, and M is the mass of 



10 



D. Sluse et al.: Size and structure of the BLR in the Einstein Cross 



microlenses. We create magnification patterns with sidelength 
of 100 r£ and pixel sizes of 0.01 r^. Tracks drawn in such a 
pattern would provide simulated microlensing lightcurves for 
one pixel-size source. To simulate lightcurves for other source 
sizes, we convolve the magnification pattern (after conversion 
on a linear flux scale) by the intensity profile of the source. As 
shown by Mortonson et al. (2005 ), the exact shape of the source 
intensity profile has little influence on the lightcurve. It instead 
depends on the characteristic size of the source. For simplicity, 
we assume a Gaussian profile. A uniform-disc profile is also 
tested for comparison. We construct magnification maps for 60 
different source sizes having characteristic scales (FWHM for 
the Gaussian and radius for the disc profile) in the range 2-1400 
pixels (0.02 r E - 14 r E ). 

4.2. Comparison with the data 

The patterns created in Sect. 14.11 are used to extract simu- 
lated microlensing lightcurves to be compared to the data. 
Two data sets are used, the OGLE data of QSO 2237+0305 
(Udalski et al. 120061) . which provide well-sampled V-band 
lightcurves of the continuum variations and our spectrophoto- 
metric lightcurves presented in Sect. [2] To study the microlens- 
ing signal, we have to get rid of intrinsic variations. This is 
done by calculating the difference lightcurves between images 
A & D. Similarly, the simulated microlensing signal is obtained 
by taking the differences between the simulated tracks for im- 
ages A & D. The parameters that characterise a simulated mi- 
crolensing track are the starting coordinates in the magnifica- 
tion patterns (xo,a, xo,d), the track position angle (9a, 0d), and 
the track length. Because the extracted track is compared to a 
lightcurve obtained over a given time range, the track length 
corresponds to the transverse velocity of the source expressed 
in r E per Julian day. As in Paper II, we assume that the velocity 
is the same for the trajectories in A & D and the track orien- 
tation 9a- 0d- To account for the roughly orthogonal shears in 
A & D (Witt & Mao 11994b , we rotate the magnification map 
of D by 90° prior to the track extraction. The difference of the 
magnitude between track A and track D should, on average, 
be equal to the one of the macro-model. Like other authors, 
we account for the uncertainty on the amount of differential 
extinction between A & D and on the macro-model flux ratio 
by allowing for a magnitude offset mo between the simulated 
lightcurves extracted for A and D. The other parameters char- 
acterising the lightcurves are the convergence and shear (k,j) 
fixed in Sect. 14. II Physical quantities (source size and velocity) 
are proportional to the mass of microlenses as M 1//2 . 

To build a representative ensemble of lightcurves that is 
compatible with the data, we followed a four-step strategy sim- 
ilar to the one described in Paper II, in Anguita et al. (|2008) 
and in Kochanek (2004): 1) we pick a set of starting values 
for the parameters defining the tracks in A & D and vary them 
to minimise the x 2 between the simulated lightcurve and the 
OGLE lightcurve; 2) we repeat step 1 to get a representative 
ensemble of 10000 tracks and good coverage of the parame- 
ter space; 3) each track estimated in step 2 is computed for 
other source sizes, and the agreement with the BLR lightcurve 
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Fig. 8. Top: Example of a good track fitting the differential 
lightcurve Aihad observed for the C IV broad line and for the 
continuum. Centre: Corresponding tracks (solid line) in the mi- 
crolensing pattern of image A (left) and D (right). The large cir- 
cle represents the C IV source size and the small one the contin- 
uum. Bottom: Track corresponding to the continuum emission 
(red solid line) and to the C IV emitting region (dashed grey 
line) for image A (left) and D (right). 



is quantified using a;f 2 -type merit function; 4) the^ 2 estimated 
in steps (2) and (3) are summed up and used to calculate a like- 
lihood distribution for any track parameter. A summary of the 
technical details regarding microlensing simulations is given in 
Appendix|B] An example of a good simulated lightcurve fitting 
the data is shown in Fig. [8] In the next subsection, we explain 
how we derive probability distributions for the quantities of in- 
terest. 

4.3. Probability distribution 

The probability of a trajectory j, defined by the set of param- 
eters p - (k, y, M, R C ,,R' S , V, ?wo,xo,a,xo,d) (where we drop the 
index j for legibility, R c s and R[, refer to the source size of the 
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continuum and of the line) given the data D is, according to 
Bayes theorem: 



P( Pj \D) = 



L(D\ P j)P(pj) 



LmpjWRcpPiVj) 



£} =1 L(D\ Pj )P(pj) 2jU L(D\ Pj )P(R c s j )P{Vj) 



,(4) 



where n is the total number of trajectories in the library, L{D\pj) 
is a likelihood estimator of the data D given the parameters pj 
and PiRl ,) and P(V7) are the priors on the continuum source 
size and velocity respectively. Only the prior on R c s and on V 
appear in Eq. |4]because we use an uniform prior on the other 
(nuisance) parameters. The prior P(R C S ■) is required because we 
have a non uniform sampling of R c s . The prior is related to the 
density of trajectories as a function of R c . given our source 
sampling in 60 different bins. If R c s . falls in the bin b of width 
lb, then we have P(R s j) °c If,. The velocity prior P{Vj) is identi- 
cal to the one used in Anguita et al. (120081 1 . It is based on previ- 
ous studies by Gil-Merino et al. (2005), which suggest that the 
net transverse velocity is lower than 0.001 r^/jd (at 90% con- 
fidence level). This corresponds to v< 660 kms 1 in the lens 
plane (~ 6000 kms~' in the source plane) for M ~ O.1M0, 
in agreement with the velocity prior used by Kochanek (2004) 
and Poindexter & Kochanek (2010a]!. 

As explained in Kochanek (2004), the use of the standard 
likelihood estimator L(D\pj) = exp(-x 2 (pj)/2) can lead to 
misleading results because it is very sensitive to the accuracy of 
the error bars and to change of the x 2 induced by possible out- 
liers. To circumvent this problem, we follow Kochanek (2004) 
and use the following likelihood estimator: 



L(D\pj) = r 



ridof-2 X 2 (Pj) 



(5) 



where F is the incomplete gamma function, and x — X ogle + 
X 2 - is the sum of the x 2 obtained when fitting the OGLE data 
and emission line lightcurves (steps 2 & 3 in Sect. I4.21 i. 

5. Inferences of the quasar structure 

In this section, we present the results of our microlensing simu- 
lations and discuss them in light of the phenomenological study 
of the emission lines in A & D. First, we discuss the domain of 
validity of our simulations and the meaning of the BLR size 
estimated with our method. Second, we derive the size of the 
continuum and of the BLR for C IV and C III]. Third, we com- 
pare our results to the radius(BLR)-luminosity relation derived 
from reverberation mapping. Fourth, we derive the black hole 
mass of QSO 2237+0305 using the virial theorem and com- 
pare it to black hole mass derived from the accretion disc size. 
Finally, we discuss the consequences of our analysis for the 
BLR structure of QSO 2237+0305. 

5.1. Caveat 

In microlensing simulations, we represent the projected BLR 
with a face-on disc having a Gaussian or a uniform intensity. 
The exact intensity profile should not affect our estimate of the 
size (Mortonson et al. 2005 ), but the sizes retrieved are only 
robust if the BLR actually has a morphology that is disc-like 



when projected on the plane of the sky. This corresponds to the 
case of a spherically symmetric BLR geometry or of an axi- 
symmetric geometry with the axis roughly pointing towards 
us. Since the microlensing signal remains nearly unchanged for 
sources with ellipticity e = 1— q < 0.3 (Fluke & Webster 1999 
Congdon et al. 2007), our simulations should also be valid for 
disc-like sources inclined by up to <p ~ 45° with respect to the 
line of sight. This means that inclination may affect our size 
estimates by a factor cos <p > 1 / V2. Although the exact ge- 
ometry of the BLR is unknown, there are several lines of evi- 
dence that QSO 2237+0305 is seen close to face-on. Firstly, 
QSO 2237+0305 is a Type 1 AGN with a flat radio spec- 
trum (a = -0.18, Falco et al. 119961) . indicating that the source 
is seen nearly face-on. Second, the analysis of the microlens- 
induced continuum flux variations by Poindexter & Kochanek 
(2010b) show that the inclination of the accretion disc is <45° 
at 68.3% confidence level. Third, the CIV line shows only a 
small blueshift that can be interpreted as evidence that the BLR 
is seen face-on (Richards et al. 120021 ). 

If the BLR has a disc-like morphology, the projected emis- 
sion should take place in a ring region (with an inner ra- 
dius Rj„ and outer radius R oul ) rather than in a uniform disc as 
assumed here. The effect of such a 'hole' in the BLR emis- 
sion has been studied by Fluke & Webster ( 1 19991 ). who show 
that the microlensing lightcurve would be significantly affected 
only for /?,„ > 0.5 R out . Since the BLR is commonly assumed 
to have R out > 10R in (e.g. Murray & Chiang [19951 Collin & 
I lure 2001 . Borguet & Hutsemekers l2010l l. we do not consider 
this geometrical issue. 

5.2. Absolute and relative size of the continuum and 
the BLR 

The V-band OGLE lightcurve is associated with emission aris- 
ing mostly from the quasar continuum. The probability distri- 
bution of the continuum source size R c s derived from our mi- 
crolensing study of that lightcurve (Sect. |4]i is displayed in 
Fig. [9] In this figure, we also show the average velocity of 
the tracks in each source size bin on a colour-coded scale. We 
see that the source size and the velocity are strongly corre- 
lated as expected from the source size/velocity degeneracy (e.g. 
Kochanek 2004 Paper II). When we consider a uniform prior 
on the velocity (not shown), we find that by imposing mo, the 
average difference of magnitude between A & D lightcurves, 
and the macro-model expectation to be equal, we exclude very 
high velocities (v > 0.005r£/jd). For the velocity range allowed 
by the Gil-Merino et al. (2005 ) velocity prior, the same distri- 
butions of mo are found regardless of what velocity bin is con- 
sidered. The median half-light radius of the continuum we de- 
rive from our simulations is R c s = 8.8 x 10 15 cm (2.5 x 10 cm 

< R c s < 25 .4x 10 15 cm at 68.3% confidence level, M = 0.3M o ) 
for a Gaussian profile, and R c s = 12.3 x 10 15 cm (3.2 x 10 15 cm 

< R c s < 36.8 x 10 15 cm at 68.3% confidence level, M = 0.3M o ) 
for a uniform disc. Table [4] compares our measurements with 



1 In the case of emission coming from a spherical shell, the pro- 
jected profile may also look ring-like, but emission should also arise 
from the central part of the disc. 



12 



D. Sluse et al.: Size and structure of the BLR in the Einstein Cross 



Table 4. Comparison of V-band continuum source size of QSO 2237+0305 (Sect. lBT2"l i. The reference is presented in Col. #1, the 
source profile used in Col. #2, and the published size estimate in Col. #4 (where R= radius of a uniform disc, Ry= scale radius 
of a Shakura-Sunyaev accretion disc, cr v and FWHM are resp. the standard deviation and the full width at half maximum of a 
Gaussian). These quantities have been converted to half-light radius R\/2 to allow comparison between profiles. Col. #6 & #7 
provide the average mass of microlenses ((M)) used in the microlensing simulations and information regarding the velocity prior. 



Reference 



Profile Estimator Size (10 15 cm) R [/2 (10 15 cm) 



(M) 



17 - 6 -i25 (FWHM) 
17.4*^ (R) 



Velocity prior 



This paper 
This paper 

Poindexter et al. J2010bt 
Anguita et al. ( 12008b 
Anguita et al. J20081 
Paper II 
Paper II 

Kochanek l l2004l > 
Kochanek l l2004l > 



Gauss 
Uniform disc 
Shakura-Sunyaev 

Gauss 

Gauss 

Gauss 

Gauss 
Shakura-Sunyaev 

Gauss 



Median 
Median 
Median 
Most likely 
Most likely 
Most likely 
Most likely 
Median 
Median 



5.o_ 23 

^+3.4 



+3 S (Rv) 



4.6*^ (TV) 
1.3*jg (cry) 
AQ%\ (FWHM) 
9.21f | (FWHM) 
4.1 + ^ (,R V ) 

5.1!^ (cry) 



°-°-6.3 
(9 t+24.5 
iz..j_ 9 ! 

14 I 92 

5 4 +4 -° 

1.5 

20+"" 5 
5 9 + " 

J -3.2 



•+0.3 
0.8 



^M/0.3M Q Gil-Merino et al. < t2005b 

^M/0.3M Q Gil-Merino et al. d2005b 

0.3 Mq Poindexter et al. d2010at 

■sjM/0.1M Q Uniform 

y/M/0AMQ Gil-Merino et al. d2005l 

^M/0AM Q Uniform 

0. 1 M Q Similar to Kochanek ( [20041 

[0. 1 , 1 ] M Q Sect. 2. 3 of Kochanek ( [20041 

[0.1,l]M o Sect. 2.3 of Kochanek ([20041 



the results of Kochanek (2004, KOC04), Anguita et al. (2008), 
Paper II, and Poindexter & Kochanek (I2010bl >. The half-light 
radius (R 1/2) allows the comparison of sizes derived with var- 
ious intensity profiles. We see that our results are compatible 
with all other works. 

The probability distributions of the half-light radius mea- 
sured for the C IV and C III] emission lines (using a Gaussian 
source light profile) are displayed in Figure [10] Reasonably 
good fits are obtained for all the line components. However, 
the quality of the best fits tends to be poorer for the broad- 
est line component (typical reduced x 2 ~ 2). The poorer fits 
were obtained for the very broad component of the C III] line 
that had best reduced x 1 ~ 4. This suggests that our error 
bars on the flux ratio are underestimated (by a factor ~ V2) 
when measured in the broadest line components. This is ex- 
pected since the flux in these components of the emission lines 
are more prone to measurement errors in the MCD decompo- 
sition. Our results are, however, robust in the sense that they 
remain unchanged when the error bars are increased by this 
factor. Median half-light radius and 68.3% confidence intervals 
are provided in Table[5]for both Gaussian and uniform-disc in- 
tensity profiles of the BLR. These sizes are independent of the 
intensity profile used for the continuum emission. These results 
confirm that the most compact region of the quasar is also the 
one emitting the broadest components of the lines. 

While the absolute source sizes derived above depend on 
the mass of the microlenses, the relative size of the contin- 
uum and the BLR are independent of the mass of microlenses. 
Table [6] displays the ratio between the emission line and the 
continuum half-light radius size R'jR'i (assuming the same 
source profile for the continuum and BLR) at 68.3% confidence 
for the different components of the BLR identified in Sect. 12.11 
We can see that the region emitting the broadest components of 
the C IV and of the C III] lines (i.e. C IV BC and C III] VBC) 
are about 4 times larger than the region emitting the V-band 
continuum. In contrast, the narrowest component of the emis- 
sion lines arises from a region that is typically 25 times larger 
than the continuum. 




B$ ( <M/0.3M e > 1/2 cm) 

Fig. 9. Probability distribution for the half-light radius of con- 
tinuum emission in QSO 2237+0305, colour-coded as a func- 
tion of the average (lens plane) velocity of the tracks in a given 
bin. 



We notice that the radius of the extreme-UV continuum 
ionising the carbon lines (typically A < 500 A) is expected to be 
much smaller than the V-band continuum, with a dependence 
of the size on wavelength as R oc /I? (Paper II and ref. therein). 
For a standard accretion disc model, ^ = 4/3, implying that the 
UV ionising continuum is more than 7 times smaller than the 
OGLE-V band continuum. 

Finally, we note that the results of Tables [5] & [6] are not 
sensitive to possible outliers in our BLR lightcurves. We es- 
sentially retrieve the same results when the points inconsistent 
with the OGLE lightcurve (i.e. those which are marked with 
a 'x' in Fig. [TJ are excluded. Similar results, but affected by 
larger error bars, are also obtained if we do not consider the 
2007 data. 
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i?,,(<Af/0.3Af B > 1 / 2 cm) 

Fig. 10. Probability distribution for the size of the C IV {top) 
and C III] (bottom) broad emission lines. Measurements for the 
region emitting the CIV & CHI] full emission line profiles 
(thick black lines, TOT) and for their components NC, BC, and 
VBC (in order of increasing FWHM, cf. Fig. [2j, are presented. 



5.3. C IV radius-luminosity relationship 

In this section, we compare our microlensing-based estimate 
of the C IV emission line (Sect. 15.2b with the radius(BLR)- 
luminosity relationship (hereafter Rblr - L) derived with the 
reverberation mapping technique (Sect. II. 2\ . The Rblr - L 
relationship for CIV has only been obtained for six objects 
(Peterson et al. 120051 Kaspi et al. 2007). The luminosity used 
in this relationship is commonly L / |(1350A). Unfortunately, 
the latter is out of our spectral range. Instead we use 
,11^(1450 A)=10 45 53 erg/s (Assef et al. l20T0l l. which has been 
shown to be equivalent to 1^(1350 A) (Fig. 4 of Vestergaard 
& Peterson 2006). This luminosity estimate, directly mea- 
sured on our spectra, is corrected for macro-magnification and 
line-of-sight Galactic extinction. The macro-magnification is 
correct within a factor 2 (Witt & Mao 1994). Because the 
lensed images are located behind the bulge of the galaxy, it 
is reasonable to assume that the absolute extinction due to 
the lens is < 1.0 mag. The finite slit width of 0.6" and see- 
ing ~ 0.6" may lead to underestimating the total flux by a 



Table 5. Half-light radius of the line emitting region R' s (in 
10 15 jM/0.3 Mq cm) for a Gaussian (G) and a uniform disc 
(D) intensity profile. Measurements for the C IV & C III] full 
emission line profiles (TOT) and for their individual compo- 
nents NC, BC and VBC (in order of increasing FWHM, cf 
Fig. [2]), are presented. 



Line 


R's (G) 


68.3% CI 


K (D) 


68.3% CI 


CIV BC 


42.9 


[14.1,124.1] 


56.6 


[17.6,168.6] 


CIV NC 


> 82.7 




> 103.5 




CIV TOT 


170.7 


[47.7,457.6] 


276.6 


[72.0,676.0] 


C III] BC 


150.1 


[51.7,430.5] 


209.6 


[69.0,627.7] 


C III] NC 


> 97.9 




> 112.2 




CIII] VBC 


53.6 


[21.2,126.9] 


41.7 


[17.4,103.0] 


C III] TOT 


126.4 


[36.9,399.5] 


163.5 


[43.8,570.6] 



Table 6. Ratio of size R'jR'i of the BLR to continuum re- 
gion for Gaussian (G) and uniform disc (D) source profiles. 
Measurements for the region emitting the C IV & C III] full 
line profiles and for their components NC, BC and VBC (in 
order of increasing FWHM), as shown in Fig. [2] are presented. 



Line 


R' s /R c s (G) 


68.3% CI 


R\IR C s (°) 


68.3% CI 


CIV BC 


3.9 


[1.6,13.4] 


3.0 


[1.1,12.0] 


CIV NC 


27.7 


[9.0,89.3] 


22.2 


[7.1,73.8] 


C IV TOT 


18.9 


[4.6,72.3] 


17.5 


[4.5,61.2] 


C III] BC 


16.9 


[4.9,63.1] 


15.2 


[4.0,58.9] 


C III] NC 


29.4 


[10.5,103.6] 


25.6 


[8.1,82.5] 


CIII] VBC 


5.0 


[2.4,14.3] 


2.3 


[0.9,7.9] 


C III] TOT 


13.8 


[0.7,53.8] 


12.6 


[3.4,49.4] 



factor up to 2. This leads to an uncertainty factor of 4 on 
/1L /1 (1450A). The CIV size has been derived to be R < ; IV = 
x 



1 7 +2 ' 9 

1,7 -1.2 



10 17 VM/0.3M o cm = 65.6+^ y/M/0.3M e light- 
days (Table [5]). Figure QT] shows the /?blr-L diagram that in- 
cludes our measurement for QSO 2237+0305 and the analyt- 
ical relation derived by Kaspi et al. (120071 ) for their so-called 
'FITEXY' fit. Although the reverberation-mapping size may 
differ from the half-light radius by a factor up to a few depend- 
ing on the geometry and density of the BLR gas, we find re- 
markably good agreement with the R-L" relationship. 

5.4. Black hole mass estimates 

The BLR size we derive from microlensing can be used to es- 
timate a black hole mass using the virial theorem: 



m bh = 1- RblrFWH m 1 blr , 

Lr 



(6) 



where f-l is a geometric factor, which encodes assump- 
tions regarding the BLR geometry (e.g. Vestergaard & 
Peterson 2006, Collin et al. 2006), calibrated from the results 
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A L A (1350 A) [erg s -1 ] 

Fig. 11. Rciv-L diagram combining reverberation mapping 
measurement of Rqiv C x ' symbol, Peterson et al. 2005 
Kaspi et al. 2007} and our microlensing-based size for 
QSO 2237+0305. The solid line is the analytical fit to the re- 
verberation mapping data from Kaspi et al. d2007l l. which has 
a slope a = 0.52. 



3960± 180km s- 1 
2.0 + ?1 x 10 s Mq. 



of OnkerHet al. d2004l l. Using FWHM C iv 
(Assef et al. 2010), we derive M BH = 5 
Consistent measurements, but generally affected by larger er- 
ror bars, are derived using the other broad components of the 
emission lines. Hereafter we refer to this black hole mass esti- 
mate as Mg^f. As we argued earlier, inclination with respect to 
the line of sight may reduce Rblr by cos </> > 1/ y2. 

Our estimate of M^ may be uncertain by a factor of a 
few because of the variation of / from object to object. A 
value of / higher than our fiducial value / = 1.0 may be 
expected for axisymmetric geometry of the BLR seen nearly 
face-on. The study of several plausible theoretical models of 
the BLR allowed Collin et al. (120061) to show that inclina- 
tion of axisymmetric BLR models may noticeably increase 
/. Observationally, Decarlj^ et al. (120081 ) find / in the range 
[0.5,10], with higher values (/ > 4) only for objects with 
FWHM < 4000 km s" 1 . 

The use of the virial theorem to derive Mbh based on the 
CIV line is still debated. For this line, non gravitational ef- 
fects, such as obscuration and radiation pressure, are proba- 
bly affecting the CIV profile (Baskin & Laor 2005, Marziani 
et al. I20061 I. Obscuration, if present, should lead to under- 
estimating FWHM(C IV) and therefore to overestimating the 
Mg^f. The effect of radiation pressure might be subtler. While 
Marconi et al. (2008) argue that radiation pressure implies 
an underestimate of Mbh, recent calculations of Netzer & 
Marziani (2010) show that orbits of BLR clouds are affected by 
radiation pressure in such a way that the product Rblr FWHM 2 



2 Onken et al. derive / = 5.5 but using a 2 instead of FWHM 2 in the 
virial relation. The corresponding value of / when using the FWHM 
is/=l. 

3 Decarli et al. use a different definition of /: /thispaper = /D Ccar ii- 



(and therefore Mbh) will vary by at most 25% owing to radi- 
ation pressure. Our finding of roughly similar Mbh for differ- 
ent components of the line, emitted in regions differently af- 
fected by radiation pressure, further supports the idea that Mbh 
in QSO 2237+0305 is not strongly biased by nongravitational 
effects. 

A black hole mass estimate can also be derived from the 
size of the accretion disc. In the accretion disc theory, the 
radius which the disc temperature equals the photon en- 
ergy, kT = hc/A rest is given by the expression (Poindexter & 
Kochanek l20TUbb 



45G4 st M, 



B1I 



M 



16 7T 6 h c 2 



1/3 



(7) 



where Mbh is the black hole mass and M is the accretion rate. 
Rewriting this expression for the black hole mass, we get 



M 



disc 



0.68x 



R 



1/2 



3/2 



\0 l5 cm) \0.2 jjm 



1/2 
%1 



■1/2 



10^0,(8) 



where 0.1 770.1 = r\ - L/(Mc 2 ) is the radiative efficiency of 
the accretion disc, R1/2 = 2.44 Rj the half-light radius of a 
Shakura-Sunyaev accretion disc model, and L/Le the lumi- 
nosity in units of Eddington luminosity. Using the accretion 
disc half-light radius R c s = 8.8 x 10 15 cm (Sect. I3T21. we get 
m bh = 17 - 7x 10 s M (2.7xlO 8 M < Mf™ <8.7xlO 9 M 
at 68.3% confidence) if we assume that L — Le and 77 =0.1. 
Unfortunately, M^ is also prone to systematic errors and de- 
pends strongly on the exact temperature profile of the disc. The 
black hole mass derived with this method assumes a temper- 
ature profile T oc R~£ with = 3/4 (Shakura & Sunyaev 
119731 ). A different value of in the range [0.5,1] may lead to 
higher/lower M^ H C by a factor as large as ~ 5 (e.g. Poindexter 
& Kochanek 2010b ). The M d BH c is thus more uncertain and less 
accurate than M BH R . 

By comparing the isotropic luminosity L = 4.0 x 10 46 erg/s 
of QSO 2237+0305 (Agol et al. 120091 to the Eddington lumi- 
nosity L E = 1.38 x 10 3& Mbh/Mq erg/s, we derive a black hole 
mass M| H > 2.9 x 10 8 Mq, consistent with M^ R and M\ 



disc 



5.5. The structure of the BLR 

Based on the MCD decomposition of the spectra, we have 
shown that the C IV and C III] lines can be decomposed into 
two (resp. three) components of increasing width (Fig. [2]). The 
possibility that these line components correspond to physically 
distinct regions in the BLR has been discussed by several au- 
thors (e.g. Brotherton et al. 1994b Marziani et al. 120101 ) but 
remains debated. The variable amount of microlensing in these 
components allows us to shed new light on this scenario. 

The differential microlensing lightcurves observed for the 
Gaussian components of the emission line (Fig. [6]) are essen- 
tially flat, even when a substantial microlensing signal is ob- 
served in the continuum. They also show an increasing mi- 
crolensing effect in the broadest line components. Our simu- 
lations show that it is possible to reproduce both the contin- 
uum and BLR lightcurves by only changing the source size. 
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We clearly find that the most compact regions have the broad- 
est emission. This excludes BLR models for which the width 
of the emission line is proportional to the size of the emitting 
region. Our results are compatible with a variation in the BLR 
size with FWHM 2 as expected if the emitting regions are viri- 
alised. This arises in various geometries, i.e., for a spherically 
isotropic BLR, a Keplerian disc, or a disc+wind model similar 
to one proposed by Murray & Chiang (11995111997b . 

The other techniques we used to analyse the spectra pro- 
vide additional evidence that the BLR of CIV and CIII] is 
composed of several regions. The "narrow-band" technique 
(Sect. 12.31 1 reveals that Am^D is closer to the macro model 
value (i.e., Am^D ~ -0.1 ± 0.1) in the core of the emission 
lines than in the wings. But even in the central 200km s _1 part 
of the line, we find Athad — 0.6 mag, indicating that there is 
microlensed flux in the core of the line. The MmD technique 
(Sect. 12.2b provides more clues to the origin of this flux. The 
shape of Fm^ for the carbon lines smoothly increases from the 
base to the line centre, and is compatible with a single-peak 
profile. The lack of "accretion disc-like" profile with a flatter 
or depleted core (see e.g. Bon et al. 2009) does not exclude 
a disc-like geometry of the BLR because disc emission can 
also produce single-peaked profiles (Robinson 1995 Murray 
& Chiang fT997l Down & al. [20T0t . The MmD supports the use 
of a Gaussian to empirically decompose the emission lines with 
the MCD technique and a composite model of the BLR formed 
by at least 2 different regions. 

Clues to the physical conditions occurring in these re- 
gions should come from the flux ratios measured between 
lines. A detailed analysis, however, needs very high signal- 
to-noise measurements (based e.g. on an average spectrum of 
QSO 2237+0305), proper account of Galactic reddening, and 
detailed photo-ionisation models. This is beyond the scope of 
this paper and will be devoted to a future publication. We 
do, however, mention several features already visible in our 
spectrum that should provide interesting constraints on photo- 
ionisation models. First, there is no broad Si III] emission. The 
ratio Si III]/C III] is very sensitive to the electronic density 
(Keenan et al. \TW!\ Aoki & Yoshida [19991 . The absence of 
Si III] in the broad component of the emission indicates that the 
electronic density n e in the region emitting the broad compo- 
nent of C III] is lesser than in the region emitting the narrow 
portion of the lines. Second, there is no clear narrow Fe II+III 
emission, suggesting that the latter arises from the inner part 
of the BLR or from a very compact region. This contrasts 
with microlensing observed in the lensed quasar RXJ 1131- 
123 1 (Sluse et al. 2007) where narrow and broad emission were 
found for the 'optical' Fe II+III emission (A e [3000, 5000] A), 
since the UV Fe II is poorly covered by the spectra in that sys- 
tem. Third, the ratio C III]/C IV is found to be nearly the same 
in the broad and narrow components of the line (Fig. |4j. This 
advocates for a similar ionisation parameter U in these two re- 
gions (Mushotzky & Ferland |1984| l, which is a priori surpris- 
ing since U should vary with the distance. In their decompo- 
sition of the emission lines based on the principal component 
analysis technique, Brotherton et al. (1994b) faced a similar 
problem and found that CIII]/CIV measured for the narrow 
part of the line (their so-called intermediate line region) was too 



large by an order of magnitude compared to photo-ionisation 
models. We speculate that metallicity and optically thick/thin 
nature of the BLR gas are likely to play important roles in re- 
producing these ratios (Snedden & Gaskell 1999). 

Small differences between the C IV and C III] line profiles 
are visible in the total line profile (Fig. HJ. The C III] profile 
peaks nearly at the systemic redshift and is rather symmetric, 
while the CIV line shows a low blueshift and an excess of 
emission in the blue (i.e. blue asymmetry), in agreement with 
previous studies (Brotherton et al. |1994ab . Figure |4] shows that 
these differences are visible in the narrow component of the 
line. On the other hand, the very broad component of these 
lines is slightly asymmetric for both lines, showing again a 
blue excess but a more greater one for C IV. Because the in- 
trinsic absorber is superimposed on the C IV profile, it is hard 
to assess whether the broadest component of the profile is in 
fact blueshifted as is the narrow one. Finally, we mention that 
the core of the carbon lines varies less than the base of the 
line, in agreement with the finding of Wilhite (2006 ) based on 
a large sample of SDSS quasars. This finding is consistent with 
the above arguments showing that the broadest part of the line 
arises from a region closer to the continuum. 

6. Conclusions 

For the first time, we have derived the size of the broad line 
region in QSO 2237+0305 based on spectrophotometric mon- 
itoring data and consisting of 39 different epochs obtained be- 
tween Oct 2004 and Dec 2007. To reach this goal, we mea- 
sured differential lightcurves between images A & D for the 
C IV and C III] broad emission lines and compared them to mi- 
crolensing simulations. This led to determining the half-light 
radius of the CIV emitting region: Rciv ~ 66 lt-days (19 lt- 
days < Rciv < 176 lt-days at 68.3% confidence), in very good 
agreement with the R blr-L relation derived from reverberation 
mapping (Kaspi et al. 2007). The size we derived for C III] is 
R cm] ~ 49 lt-days (14 lt-days < # C iii] < 154 lt-days at 68.3% 
confidence), compatible with the size estimated for the region 
emitting C IV. 

Thanks to the variable amount of microlensing observed 
within a given emission line, we can also derive information on 
the structure of the broad line region. Differential lightcurves 
obtained for various velocity slices of the C IV and C III] lines 
show that the wings of the lines are more microlensed than the 
core, indicating that the former arise in a more compact re- 
gion. This finding is confirmed by two other techniques. The 
first one demonstrates that a broad and single-peaked fraction 
of the emission line is microlensed, while a narrower frac- 
tion is unaffected by microlensing. This technique also sug- 
gests a slightly different structure for the C IV and C III] emis- 
sion regions, with the narrow C IV emission blueshifted with 
respect to the systemic redshift. The second technique as- 
sumes that the emission lines can be decomposed into a sum 
of Gaussian components. The CIV is separated into a narrow 
(FWHM ~ 2600^8-') and a broad (FWHM ~ 6300 kms" 1 ) 
component. In order to reproduce the small differences be- 
tween the C IV and C III] lines, we need three Gaussian pro- 
files (FWHM ~ 1550, 3400, 8550km s" 1 ) to reproduce the lat- 
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ter. The lightcurves derived for these components are essen- 
tially flat and show that microlensing is more important when 
the FWHM increases. Although the individual Gaussian line 
components do not necessarily isolate individual emission re- 
gions, we find that these lightcurves are compatible with the 
continuum lightcurves, provided only the size of the emis- 
sion region is modified. This allows us to derive a half-light 
radius for the regions emitting these components, as well as 
for their size relative to the continuum. The radii are con- 
sistent with the virial hypothesis and a radius that varies as 
FWHM 2 . Using the virial theorem, we derived a black hole 
mass M BH ~2.0xl0 8 M o (0.5xl0 8 M o < M BH < 5.4xlO 8 M 
at 68.3% confidence). 

Our analysis supports the findings by other authors 
(Brotherton et al. U994al Marziani et al. 1201 Oi l that the regions 
emitting the C IV and C III] lines are composed of at least two 
spatially distinct components, one emitting the narrow core of 
the line and another, more compact, emitting a broadest com- 
ponent. The broad (resp. narrow) component of the C IV and 
C III] lines do not have exactly the same profiles. The flux ratio 
C III]/C IV is very similar when measured in the broad and in 
the narrow components of the lines. This suggests that the ioni- 
sation parameter U is nearly the same in the two regions, a sur- 
prising result since the narrow components are found to arise 
in regions that are several times larger than the broad compo- 
nents. Other lines observed in our spectra seem to arise from a 
least two components: MgIU2798, A1IU1671, HeIU1640, 
A1IIU1857. The situation is different for Si III] ^1892, which 
does not show emission from a broad component, which is in- 
dicative of a smaller electronic density n e in the region emitting 
the broadest part of the emission lines. On the other hand, we 
detect broad microlensed Fe II+III but no "extended" emission. 
This suggests that Fej/y is produced in the inner part of the 
BLR or in a very compact region. Obtaining spectrophotomet- 
ric monitoring data in the near-infrared where Balmer lines are 
detectable would be very useful for constraining photoionisa- 
tion models and comparing the microlensing signal in Balmer 
(e.g. HjS), high ionisation (e.g. C IV), and Fe op , lines. 

We demonstrated that the spectrophotometric monitoring 
of microlensing in a lensed quasar is a powerful technique for 
probing the inner regions of quasars, measuring the size of the 
broad line region, and infering its structure. More work is still 
needed to take full advantage of the method, but our results are 
very promising. Several improvements are possible to increase 
the accuracy of our size measurements and better characterise 
the BLR structure. First, microlensing simulations reproduc- 
ing the signal in more than 2 lensed images should allow one 
to narrow the final probability distribution on the source size. 
Second, the implementation of spectral fitting using Markov- 
Chain Monte-Carlo should allow a more appropriate estimate 
of the error bars and more advanced modelling of the individ- 
ual spectral components. Third, a fully coherent scheme should 
be developed to consistently model the emission line shape and 
the corresponding source intensity profile used in the simula- 
tion. These improvement will be the subject of a future work. 
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Appendix A: Averaged emission line profile 

We explain here how we construct an averaged emission line 
profile and how we use it to choose a reasonable multi- 
Gaussian decomposition of the broad emission lines. Following 
the prescriptions of the multi-component decomposition de- 
scribed in Sec. 12.11 we are able to subtract the quasar pseudo- 
continuum emission (including the Fe II+III emission) from the 
quasar spectra. This procedure leads to a spectrum containing 
only the broad emission lines. Then, we normalize their flux 
and compare the line profiles averaged over periods P 1 , P2, and 
P3. This procedure reveals that the profile of the C III] emission 
line basically remains unchanged over the whole monitoring 
period. In contrast, there is evidence of a small change in the 
C IV line profile between 2005 and 2007. The high signal- to- 
noise emission line spectra of C III] and C IV are then fitted 
with a sum of 3 Gaussian-line profiles. For C III], we find that 
the best fit is obtained with a narrow emission component with 
FWHM = 1545 km s~', an intermediate emission component 
FWHM = 3400 km s~', and a broad component with FWHM 
= 8548 km s~'. For CIV, the emission line is well fitted with 
only two emission components with FWHM ~ 2570 km s _I and 
FWHM ~ 6150km s _1 . The absorption component is also fit- 
ted with an FWHM ~ 1250km s 1 for PI and P2 and FWHM- 
900 km s 1 for P3. The difference in width of the absorption 
comes from the increased resolution of the spectra obtained 
during period P3. Higher resolution spectra (Hintzen et al. 
[19901 Yee & De Robertis [19911) show that the C IV absorp- 
tion system is composed of at least two different absorption 
clouds associated with the quasar. Narrow absorption systems 
are also present in front of C III] emission (Motta et al. 2004 ), 
but they do not appear as significant features in our spectra and 
are therefore not part of our modelling. 

If we now compare the decomposition of the line profile 
in images A & D, we find that the width of the individual 
Gaussian profiles are similar but not their relative intensities, 
as already suggested by the analysis in Paper I. During our fit- 
ting procedure, we assume that the FWHM of the Gaussian 
components are identical for A & D, except for C IV(BC) and 
C IV(AC). The widths are also fixed to the fiducial values de- 
rived here. 



ference lightcurve Am' k (p) is then compared to the data D = 
Am* by measuring the goodness of fit with ax 2 statistic, 

x\p) = Y\ — > (B - 1} 

where <7> are the uncertainties of the data and « b s is the num- 
ber of data points. For OGLE data, we use those obtained be- 
tween modified julian day (i.e. jd-2450000) 3 100 and 4500, and 
we bin the data points obtained the same night, which leads to 
"obs = 181. The error at is the quadratic sum of the photomet- 
ric error o~ ogle provided by OGLE, of the standard deviation 
a"bin between binned points, and of the systematic underestima- 
tion of the error (T sys of the OGLE data. The error <x svs is esti- 
mated using a polynomial fitting of the difference lightcurves 
(see Paper II) and is estimated to be cr sys = 0.05 mag. The 
data points corresponding to the BLR lightcurve have been de- 
scribed in Sect. 13.11 

The brute force approach, which consists in randomly pick- 
ing track parameters and estimating a x 2 is very inefficient in 
getting even a few tracks reproducing the data unless a large 
amount of computing time is available. To overcome this prob- 
lem, we follow the procedure explained in Paper II and, for 
each set of random parameters we have generated, we optimise 
the 6 parameters (mo,X\,XD,ff) with a x 2 based minimisation 
algorithm using a Levenberg-Marquardt least-square routine. 
Contrary to what was done in Paper II, we do not optimise 
the track velocity to avoid possible bias on the velocity dis- 
tribution induced by the minimisation routine. In addition, we 
windowed the magnification map in such a way that the starting 
point of a track is always more then one track length away from 
the border of the map. Our final library of tracks fitting OGLE 
data contains 6xl0 5 tracks (10000 track per source size) with a 
X 2 /n obs < 8.5. 

The only parameter that is varied to compare the tracks of 
our library to the BLR lightcurves is the source size R s . As ex- 
plained in Sect. 14.21 we extract, for each track of our library 
(step 2), the same track in the 59 other magnification maps cor- 
responding to the 59 source sizes. We then search for the best 
source size using a minimisation routine and calculate tracks 
for arbitrary source sizes (in the range 0.01-14 rg) by means of 
a quadratic spline interpolation between different source sizes. 



Appendix B: Microlensing fit 

We provide here the technical details regarding the compari- 
son of the simulated tracks with data. The following source 
sizes are considered: 2, 4, 8, 12, 16, 20, 24, 28, 32, 36, 40, 
44, 48, 52, 56, 60, 70, 80, 90, 100, 110, 120, 130, 140, 150, 
160, 170, 180, 190, 200, 240, 280, 320, 360, 400, 440, 480, 
520, 560, 600, 640, 680, 720, 760, 800, 840, 880, 920, 960, 
1000, 1040, 1080, 1120, 1160, 1200, 1240, 1280, 1320, 1360, 
and 1400 pixels. For each of the 60 source sizes R s , we sim- 
ulate 10000 pairs of lightcurves for both images A & D by 
tracing source trajectories across the magnification patterns. 
Bilinear interpolation is used for the track extraction. Each 
track is characterised by the set of parameters described in 
Sect. 14.21 p = (R s , V, mo, 0, xo,a,xo,d, M). Each simulated dif- 



